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Probing the Origins of Neutrino Mass with Supernova Data 

Hooman Davoudiasl* and Patrick Hubert 
Department of Physics, University of Wisconsin, Madison, WI 53706, USA 

We study type II supernova signatures of neutrino mass generation via symmetry breaking at a 
scale in the range from keV to MeV. The scalar responsible for symmetry breaking can be thermalized 
in the supernova core and restore the symmetry. The neutrinos from scalar decays have about half 
the average energy of thermal neutrinos and are Bose-Einstein distributed. We find that, even 
without a detailed knowledge of the supernova parameters, a discovery is well within reach at 
Super-Kamiokande. 
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Type II supernovae provide some of the most extreme 
environments at the present cosmological epoch. In 
these explosions, the stellar gravitational binding energy 
_EgN ~ 3 X 10^^ erg is released almost entirely in a neu- 
trino burst. After the initial collapse, the stellar core 
is characterized by nuclear densities and a temperature 
TsN — 30 MeV. The cooling of the core occurs over a 
period of seconds, through an explosive release of Esn 
by neutrino emission. 

Given such conditions, it is interesting to ask whether 
there may be new particles at a scale below Tsn that 
could be relevant in the first few seconds of the super- 
nova explosion. In fact, supernova cooling has provided 
stringent bounds on many models containing light par- 
ticles, like the axion However, such models do not 
generically yield distinct supernova signatures. 

In this Letter, we study possible supernova signatures 
of neutrino mass generation from low energy symme- 
try breaking. The basic idea is that, given a typical 
supernova environment with temperature TgN, a light 
scalar, whose vacuum expectation value (vEv) generates 
rrii, 7^ 0, may be in thermal equilibrium during the explo- 
sion. We show that this can typically be the case, pro- 
viding robust supernova signatures of this physics, which 
may otherwise be hard to access. 

Our study is based on recent low-cutoff-scale models 
for neutrino mass generation proposed in Ref. Q. In 
contrast to the seesaw mechanism with scales of order 
lO^-* GeV, these models can be cut off at 10 TeV < A < 
1000 TeV, avoiding extrapolations far above the exper- 
imental frontier at around 1 TeV. Both Majorana and 
Dirac masses can be generated in these theories. To af- 
ford a low cutoff scale, a discrete gauged (anomaly-free) 
Z| X Zg symmetry group is assumed to protect Baryon 
{B) and Lepton (L) numbers. Here, small values of 
nil, 7^ are generated when certain symmetries are spon- 
taneously broken at or below 10 MeV, by a scalar VEV^. 
In this work, we generically denote such a scalar by tp. 

A consequence of this scenario is that a Bose-Einstein 
gas of scalars with T ~ TgN will be present during the 
cooling of the supernova. Thermal scalars produced in 



the bulk are either confined within the neutrino-sphere or 
they decay into neutrinos, which are also confined. The 
ip's escaping at the surface quickly decay into neutrinos. 
The average energy of these neutrinos is half that of the 
thermal scalars. In addition, these neutrinos inherit the 
Bose-Einstein distribution of the parent scalars. These 
are distinct signatures of the presence of a light bosonic 
degree of freedom that couples to neutrinos. Detecting a 
fiux of supernova neutrinos with about half the expected 
average energy and a Bose-Einstein thermal origin will 
strongly favor the type of models we consider. In this 
case, the supernova data will provide a direct probe of 
the mechanism responsible for neutrino mass generation. 

Since the thermalized scalars are expected to have a 
self coupling A ~ 1, the corrections to their potentials 
are of order XT^ :s> ((ys)^, where T ~ Tsn- Thus, another 
consequence of these models is that the symmetry whose 
breaking gives m^, 7^ is restored in the supernova core. 
We then expect that (f) = = in the core, during 
the initial stages of the supernova explosion. Thus a sym- 
metry breaking phase transition will take place after the 
star has cooled to small enough temperatures, T ^ (ip). 
This will result in the emission of non-thermal neutrinos 
with El, ~ (if). Later, we will briefiy comment on the 
possibility of detecting such a signal. 

To demonstrate these features, we adopt a neutrino 
mass model with low scale symmetry breaking, presented 
in Ref. In this model, Majorana neutrino masses con- 
sistent with flavor oscillation data are generated. Other 
models with more degrees of freedom that lead to Dirac 
masses from low scale symmetry breaking have also been 
considered Q . To keep the discussion simple and empha- 
size the main features, we will not consider these latter 
models here. However, a similar analysis can be per- 
formed for the Dirac mass models of Ref. §. 

Without specific assumptions about the cutoff scale A, 
a gauged Z| x Z| protects B and L numbers at safe 
levels Q]. To obtain Majorana masses for neutrinos, a 
scalar ip with Z| charge -1-1 is introduced. Suppressing 
the generation indices, the dimension-6 operator 
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^ Refs.p, b| contain similar ideas, using global symmetries. 



then yields the Yukawa coupling y^pvLi^Li with y = 
{v/h-Y and V = {H) =^ 0; L and H are the Standard 
Model lepton doublet and the Higgs, respectively. 



The coupling y is constrained by Qvvf3 + (/? Q : y < 
3 X 10-^ implying that A > 30TeV. For ~ 0.1 eV, 
we then need {ip) ^ 3 keV. On the other hand, bounds 
on cosmological domain walls from broken discrete sym- 
metries require {if) < 1 MeV, which yields A ^ 10'^ TeV. 
In this work, we will then consider 

1 keV < ((ys) < 1 MeV . (2) 

The upper-bound on ((^) can easily be increased by 0(1) 
factors if the domain wall network is frustrated or if the 
reheat temperature in the Early Universe is below the 
symmetry breaking scale. 

In a typical supernova environment, the scalar ip will 
be thermally produced in equilibrium and, consequently, 
the broken symmetry responsible for rriu ^ will be re- 
stored. To see this, we first note that after the initial 
collapse of the star, the hot and dense gas of particles 
around the core of the progenitor traps neutrinos, form- 
ing a neutrino- sphere. The i/-confinement typically lasts 
over the cooling time, Atcooi ~ 10 s, of the supernova. 
Assuming that the neutrinos are in thermal equilibrium 
inside the neutrino-sphere, their number density is given 
by riv ~ TgN- rate of the reaction vlvl ^ is given 
by ~ y^Ta^n which, up to phase space factors, is of 
the same order as the rate for vlul f^^, given a ip 
4-point-coupling A 2± 1.^ From (^), we expect y > 10^^, 
with ~ 0.1 eV. This implies F^ > (10-Ss)"^ 

Thus, F^^ <^ Aicooi- Note that F^^ also sets the rate for 
v^ip and v^Lp v^pp, mediated by the i-channel 

exchange of p. The mean-free-path dip for the interac- 
tions of the (^'s produced inside the neutrino-sphere is 
then given by F~^ < Im. The size of the 

neutrino-sphere is comparable to that of a neutron-star 
and we have R^, ~ 50 km 3> di^. Given all these con- 
siderations, we conclude that the neutrinos and the (/?'s 
come to thermal equilibrium inside the neutrino-sphere, 
well before the star begins cooling. This means that the 
number density of p^s is given by n^p ^ ^ ^sn- 

To study finite temperature effects, we write down a 
Z|-invariant potential for p below the weak scale 

V{p) = ypvL^L - [^w"^ + "^^^ + '^(w^)^ + h.c. (3) 

In a typical theory, /i ~ m ~ ((/?) , and A ~ 1 . Since ex- 
perimental constraints require ?/ < 3 x 10^^, the Yukawa 
coupling is much smaller than other typical couplings of 
the p system. Hence, Yukawa contributions to thermal 
corrections of are ignored in our analysis. 

Once in equilibrium, the p field obtains a thermal mass 
~ AT^. At r ~ TsN, this correction is much larger than 
the typical negative mass squared responsible for sym- 
metry breaking: AT|j^ :i> Hence, the spontaneously 
broken symmetry is restored within the neutrino-sphere. 



There is also the process vi,vi^ ipip mediated by the coupUng 
mip^ . With a natural potential for ip, we have m ~ rriip ~ {ip). 
For T 2> (<fi), this process is then sub-dominant to the first two. 



The temperature of the neutrino-sphere will eventu- 
ally fall below the critical temperature Tc at which the 
symmetry is broken. Here, we note that due to the pres- 
ence of the term mp^ , the transition back to the broken 
phase may be first order. Quite generically, we expect 
Tc — (p). Once T < Tc, the field p will roll or tunnel to 
its vacuum value. In doing so, p will oscillate about its 
VEV and radiate its energy in neutrinos. These neutrinos 
will carry a typical energy of order ~ {p) . 

To estimate the potential signal, we first note that the 
total energy stored in the symmetric phase is given 
by (p)'^ V; V is the volume of the neutrino-sphere. 
Comparing to Esn for {p) = 1 MeV, we see that only 
about 10^^ i?sN will be emitted directly during the phase 
transition. To have a chance to detect such a tiny neu- 
trino emission, only the largest available detector with 
a very low threshold can be used. The IceCube detec- 
tor has excellent capabilities to study supernova neutri- 
nos, although its single event energy threshold is around 
100 GeV. Instead, supernova neutrinos will be detected 
only by the ice glow their charged current interactions 
produce, which leads to an increase in the noise rate of 
the photomultiplier tubes The signal in IceCube is 
thus proportional to E^, where a factor E^ is due to 
the energy dependence of the cross section and an addi- 
tional factor El, accounts for the fact that the number 
of Cerenkov photons is directly proportional to Ei,. The 
number of neutrinos emitted in the transition scales as 
(p)'^- Since E^, ~ (p), the total signal scales like (p)^- On 
the other hand, the minimum signal strength needed for 
a significant detection depends on the noise of the photo- 
tubes / and on the time interval of the neutrino emission 
At^p like 1 / yJfKUp. Using a realistic set of numbers to 
describe IceCube from Ref. j6|, we find that a minimal 
signal for the phase transition requires {p) > 5 MeV and 
At^ < IQ-'^ s. This would require that the usual upper 
bounds on (p) be somewhat relaxed and also that the 
whole neutrino-sphere undergo the phase transition at 
the same time, in a neutrino-flash. 

Next, we will present our quantitative results for the 
expected size of the signal from p-g&s decays at the 
neutrino-sphere. In thermal equilibrium, p's will track 
the neutrino density and temperature closely, ie. they are 
confined within the neutrino-sphere by the p vv reac- 
tions. Neutrinos escape from the surface of the neutrino- 
sphere. Similarly, p^s also escape, but decay into neutri- 
nos promptly, as they have a lifetime r ~ T/{Xm^)'^ ~ 
10"^ s, in the star's rest-frame. Thus, the neutrino- 
sphere will radiate Fermi-Dirac thermal neutrinos, as well 
as those from decays of the Bose-Einstein (ys-gas. 

The relative luminosity is then determined by the en- 
ergy fractions carried by neutrinos and p^s. The two- 
component gas of neutrinos and p's will have its energy 
equipartitioned into all available degrees of freedom: 3 
fermionic and 1 bosonic. Thus, the total number of effec- 
tively massless degrees of freedom is 3 x (7/8) -I- 1 = 29/8. 
Hence, a fraction Fb = 8/29 ~ 0.28 of the total neutrino 
luminosity will be emitted from a Bose-distributed gas. 
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As the most important detection channel is inverse f3- 
decay, the next question is: how large is the i^j-decay 
contribution to the flux? 

The decay of the ip particles produces mass eigenstates 
instead of flavor eigenstates, hence the branching ratio 
into a flavor a depends on the hierarchy of masses rrii 
and the relevant mixing matrix elements Uai '■ 

BrK) a^|[/„,|2m2. (4) 

i 

Here, we note that neutrino oscillation data have pro- 
vided some information on Uai 0: |t^ei|^ — 0.7 and 
C/e2|^ — 0.3 whereas for |?7e3|^ there is only an upper- 
bound of 0.05 at the 3a level. Since Br(i^Q,) is propor- 
tional to mf, we can distinguish three cases. The first 
one is normal hierarchy, ie. ms is much heavier than 
mi and 1712 ■ In this case basically only the state ma 
is produced and Br(Pe) = 1 x lUes]"^ < 0.05. In the 
case of inverted hierarchy, mi and m2 have nearly the 
same mass and are much heavier than m^, which gives 
Br(Pe) = (1/2) X |t/ei|2 + (1/2) X \Ue2\^ ^ 1/2. In the 
case of degenerate neutrinos, all three masses are approx- 
imately the same and Bi^Ve) — (1/3) x |J7eiP + (1/3) x 
\Ue2\'^ + (1/3) X |[/e3|2 ~ 1/3. The branching ratio of the 
(/?-decay determines the flavor ratio for the bosonic flux 
component: 77^ — Br(i/(j) . 

In a simple blackbody approximation to the neutrino 
emission from the proto- neutron star, there would be 
again equipartition and thus the flavor ratio for the 
Fermi-Dirac component would be = 1/3. Certainly, 
the effect in Pe would be very small for normal hierarchy 
in the absence of neutrino oscillations. More detailed 
studies of neutrino emission including various levels of 
microphysical detail show a more or less pronounced dif- 
ference from the blackbody Ansatz and deviations in rj^^ 
from 1/3 up to a factor of two seem possible 

In general, the Pe flux arriving at the detector is a 
sum of 6 different initial fluxes at the supernova 

^,^^NxY, Pe.s {FdvZ + PbvI 4't ) , (5) 

Pa 

where Pgg denotes the probability of neutrinos of flavor 
a to arrive as De at the detector. is a normalization 
factor accounting for the energy and distance of the su- 
pernova. Here, there is basically no difference between 
and Vt. Both interact with the neutron star material 
only via neutral current interaction and any differences 
in their initial spectrum would be eliminated by neutrino 
oscillation via the atmospheric angle. Thus, it is suffi- 
cient to treat and v-r as one flavor 9^: with twice the 
flux. Having now an effective two flavor problem we can 
write Pxe as 1 — Pgg. Since Fb + Fd ~ 1, Eq. (||) yields 

- A^ X {Fgg [(1 - FB)TlZ<t>Z + FbiiIc^I] 

+ (1 - P,,) [(1 - Fb)vIc^I + FbvI^I] } -(6) 

The ij^^s are determined by the solar mixing angle and 
the mass ordering. Therefore, we will discuss the different 



possibilities for the mass ordering. However, within each 
case, we will assume each t]^ to be known, since the solar 
mixing angle is well constrained. Pgg also depends on the 
solar mixing angle, the mass ordering, and additionally 
on the small angle ^13. The value of Pgg is determined 
by the adiabaticity of the neutrino evolution at the MSW 
resonance which corresponds to the larger, atmospheric 
mass splitting. The adiabaticity itself depends on the 
value of ^la; for a detailed derivation see eg. |^. 

Since there is only an upper-bound on ^13 , we will con- 
sider two extreme cases: either sin^ ^13 is much smaller 
than 10"'^ or much larger. In the first case, Pgg = 
cos^ 9i2 — 0.3, irrespective of the mass ordering; in the 
second case Pgg = cos^ 612 for normal mass hierarchy and 
Pgg = for inverted hierarchy. Thus, at least a fraction 
1 — cos^ 9i2 ~ 0.7 of the i>x flux will contribute to the ob- 
servable Ve flux. Therefore, there is always a sizable con- 
tribution to the signal by neutrinos from (/3-decay. Even 
for normal hierarchy, where nearly all (/3's decay into D^, 
since at least 70% of Dx will oscillate into D^. 

The spectral distribution for the initial fluxes in the 
Fermi-Dirac case is given by the usual expression. How- 
ever, for the neutrinos originating in (p-decays, we have to 
consider that each ip particle will decay into two neutri- 
nos, with flat energy distributions in the star's rest-frame. 
Hence, the distribution functions we use are 

and 

*£(i=;r.) = 2/;f (-^), (8) 

where the factor of 2 in (f>^ accounts for 2 neutrinos per 
decay. We implicitly assume that the temperature of the 
f-gSiS and the Fermi-Dirac gas are the same, but we al- 
low for different temperatures of i>e and Dx ■ The signature 
for the presence of a </?-gas is thus the observation of up 
to four different distributions which compose the signal. 
Two of them are the usual Fermi-Dirac distributed De 
and i^x contributions, and we have another two contribu- 
tions which result from iy9-decay. The latter two have an 
average energy which is approximately only 1/2 of that 
corresponding to the supernova temperature. 

Our simple Ansatz does not include possible changes 
in the spectral shape with respect to the pure Fermi- 
Dirac or Bose distributions. Those changes arise due 
to the energy dependent position of the neutrino-sphere 
and due to diffusion processes on the way out from the 
neutron star. Since neutrino density is relatively low 
outside the neutrino-sphere, Pauli-blocking can be ig- 
nored Therefore, one naively expects that the Bose 
and Fermi-Dirac components do not interact on their way 
out. Thus, transport or diffusion should affect both com- 
ponents in a similar fashion. To understand the effects of 
the (/3-gas on the position of the neutrino-sphere(s) and its 
energy and flavor dependence requires a detailed calcula- 
tion along the lines of Ref . S and is beyond the scope of 
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this Letter. However, it seems unlikely that such effects 
could conspire to destroy the four-component feature. 

It remains to estimate the size of the event sample 
needed to detect the bi-modal nature of the v,, spec- 
trum. To this end, we convert the flux in Eq. (|^) 
to an event rate by convolving it with the cross sec- 
tion for inverse /3-decay which is approximately given 
by cr = 9.52 x lO"'*'* cm^ MeV"^ with a Q-value of 
1.8 MeV. We can now compute the event spectrum for 
a given set of parameters Tp^, T^^, 77^, Pgg, -Fb and ask 
how well we can determine or constrain Fb from a fit to 
these simulated data. For the actual computation we will 
use the inverse /3-decay cross section from Ref. [p^ and 
assume a typical detector threshold of 7 MeV We 
use the supernova input parameters 

Tp^ = 4.8 MeV, T^^ = 4.8 MeV, 77^^ = 0.66 , (9) 

corresponding to the simulation in Ref. ||^ . With Pgg = 
cos^ 6*12, a supernova with a total energy release of 3 x 
10^^ erg at a distance of lOkpc, yields ^ 10 000 events, 
in a detector with a fiducial mass of 22.5 kt, like Super- 
Kamiokande. We simulate our data assuming that indeed 
there is a iy9-gas, ie. Fb = 0.28. 

We perform a least square fit to the simulated data in 
the usual way by constructing a x^-function and asking 
how many events are needed to exclude the absence of the 
f-gds, Fb = 0, at 3 (T confidence level. More technically, 
we want to know how large Ax^(-Fb = 0) is per event. 
The fit is constrained to Tp^,!"?, > 3.8 MeV, which is 
a conservative range The number of events needed 
for a 3 cr discovery, under various assumptions about our 
knowledge of the supernova, is shown in the table. From 
this table, we see that Super-Kamiokandc is quite capable 
of detecting the presence of the tp-gas without the need 
for a detailed model of the supernova physics. 

Establishing the Bose-Einstein origin of the flux from 
(/3-decay will require a large data sample. Here, one can 
perform an analysis similar to the one above, using a 
Fermi-Dirac distribution for the scalar gas when testing 
the difference in the fit to the simulated data. We find 
that 10^ - 10^ events are needed to reach a conclusion at 



3 cr. For this purpose, our estimate suggests, one of the 
megaton water Cerenkov detectors will be needed. 
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Inverted 


P-e-e 


COS^ 6*12 


COS^ 012 





Perfect knowledge 


538 


152 


294 


N unknown 


1107 


678 


682 


N and Tp^ unknown 


1157 


709 


682 


Tp^ and Tp^ unknown 


10502 


7397 


709 


N, Tp^ , Tp^ and ry? unknown 


13449 


7432 


709 



In this Letter, we presented supernova signatures of 
neutrino mass models where new massive scalars with 
VEv's in the range keV-MeV generate mp ^ 0. We 
showed that the scalars will typically come to thermal 
equilibrium during the initial stages of supernova cool- 
ing. The subsequent decay of the scalars yields a neu- 
trino flux with roughly half the average energy of ther- 
mally emitted neutrinos. The lower energy component 
of the neutrino flux encodes the Bose distribution of the 
parent scalars. A bi-modal energy distribution for each 
flavor is a distinct and robust signature of the models we 
consider. This effect could be readily detectable at the 
3 cr level at Super-Kamiokande, for a typical galactic su- 
pernova. Observation of these signatures may offer the 
only direct access to the physics responsible for neutrino 
mass generation. 
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